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. METHOD OF PROCESSING MARINE MAGNETIC GRADIENT DATA AMD 
EXPLORATION METHODS USING THAT DATA 

Field of the Invention 1 
This invention relates to a method of processing marine 
magnetic gradient data and also to exploration methods 
using that data. 

Background Art 

Exploration of the earth to locate oil, minerals or other 
valuable deposits, uses seismic techniques in order to 
provide information about the subsurface structure of the 
earth so that an assessment can be made as to whether a 
particular deposit may or may not be present. in 
particular, seismic techniques provide an indication of 
various subsurface structures including porous sandstone 
or fractured carbonates which may contain oil or other 
valuable deposits. 



Whilst this technique does provide an indication of 
subsurface structures, it is possible a particular 
structure which may appear to be of interest, in fact 
represents a significant hazard to drilling. For example, 
it is possible that subsurface volcanos exhibit very 
25 similar characteristics to anticlines which may contain 

oil. if a drilling operation is initiated and it is found 
that the anticline is in fact a volcano, then a 
considerable financial loss results because of the cost of 
establishing the drilling platform. This is particularly 
3 0 the case in relation to marine exploration because of the 
significantly higher drilling cost. 

in order to provide more information concerning the nature 
of subsurface structures, magnetic data of the exploration 
35 area is obtained. 
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A magnetic gradient survey enables magnetic data relating 
to the survey region to be obtained and this magnetic data 
can be used to provide information relating to the nature 
of subsurface structures. if the magnetic survey is 
overlayed with the seismic survey, then structures which 
appear to be of interest from the point of view of the 
seismic survey, can be further considered in the light of 
the magnetic data so that a clearer indication can be 
formed as to whether the structure is a structure which 
may contain a deposit of interest such as an oil deposit, 
or whether the structure exhibits magnetic phenomena, such 
would be the case with a volcano. Thus, the location of 
drilling platforms can be decided with more precision to 
avoid subsurface structures which could be hazardous from 
the point of view of a drilling operation. However, 
conventional techniques for processing magnetic gradient 
data contain considerable distortion and anomalies due to 
unwanted magnetic effects including the ship bias referred 
to above. 

The conventional method of obtaining that data in relation 
to onshore exploration is to tow magnetometers behind an 
aircraft to obtain magnetic data in relation to the survey 
area. In the case of marine exploration, the 
magnetometers are towed behind a ship. 

The usefulness of marine magnetic data which is obtained 
in this manner has been limited by the quality of the data 
which is obtained. One of the major problems with marine 
magnetic data acquisition is interference or so-called 
ship bias which is created by the magnetic field induced 
by the ship which tows the magnetometers. 

In conventional marine magnetic data acquisition, two 
magnetic field sensors, which are generally referred to as 
fish, are towed behind a ship. The magnetometers are 
connected to a tow line and the magnetometer which is 

H ! xu.l fl aXK«pxSpeciNBHP Prov - netted ot Processl^ H«i 00 Miotic Cr.dl.nt O-t-.doc J/12,01 



- 4 - 



10 



15 



closest to the ship is towed at a distance of some 300-600 
metres behind the ship to avoid the effect of ship induced 
magnetic field. Furthermore, the magnetometers are 
separated by a distance of in excess of 100 metres. The 
reason for the length of the tow line and the separation 
of the magnetometers is to reduce ship bias and therefore 
provide data which is relatively free of that bias. 

However, because the towing line is so long and the 
distance of separation so great, the degree of drift of 
the magnetometers in the ocean as the magnetometers are 
towed behind the ship is considerable. Furthermore, the 
processing techniques used to acquire magnetic gradient 
data assume that after a given time period, the trailing 
magnetometer will be towed to a position which coincides 
with a previous position of the magnetometer which is 
closest to the ship. This assumption is made during the 
processing of the magnetic data. However, because the tow 
lxne is extremely long and the distance of separation 
between the magnetometers is considerable, drift of the 
magnetometers due to sea currents and the like means it is 
unlikely that the trailing magnetometer will actually 
occupy the same position as the first magnetometer after 
that given time period. 

Furthermore, in the acquisition of the magnetic data, the 
ship is required to travel along predetermined survey 
lines and because of the drifting of the sensors, the 
sensors do not actually travel along the survey line. m 
the southern hemisphere, if the magnetometers drift to the 
north of the survey line, the recorded gradient data 
between the two magnetometers will have larger distortion 
than if drifting towards the south of the survey line. 

35 Summary of the Invention 

The object of the invention is to provide an improved 
processing method which accounts for ship-induced bias and 
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which therefore enables the magnetometers to be towed by a 
tow line which is much shorter than the conventional tow 
line and also be located closer together to reduce the 
amount of drift and also increase the likelihood that the 
trailing magnetometer will occupy the same position which 
was previously occupied by the leading magnetometer after 
a given time period. 



10 



The invention may be said to reside in a method of 
processing marine magnetic data obtained by towing first 
and second spaced apart sensors behind a ship, the method 
including: 

obtaining raw magnetic gradient data from the 

sensors ; 

15 determining the trend of the gradient of the ship 

bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
20 a data output. 

Because the method determines the trend of the gradient of 
the ship bias, a more accurate estimation of the actual 
gradient bias introduced by the ship, having regard to any 
drifting of the sensors, is taken into account. By 
subtracting the trend of the gradient bias from the raw 
gradient data, the bias introduced by the ship can be 
removed from the data so that more accurate data is 
obtained. Because the present method eliminates the 
effect of the ship bias, the sensors can be towed much 
closer to the ship and at much closer separation 
distances, thereby reducing the amount of drift of the 
sensors and therefore resulting in the sensors being more 
likely to follow the actual survey line upon which the 
ship travels and also for the rear sensor to occupy the 
same position at a later time which was previously 
occupied by the front sensor. Thus, the quality of the 
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output data is greatly improved and therefore provides 
data which gives a much more accurate indication of the 
subsurface magnetic characteristics of the survey area, 
which can be then used in combination with seismic data to 
5 make an assessment as to the viability of further 
exploration or a drilling program. 

Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
10 magnetic gradient data obtained by the sensors, and 

determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

The estimate of the gradient of the ship bias may in fact 
15 simply be the raw magnetic gradient data which is obtained 
from the sensors which may include data relating to the 
gradient of the ship bias as well as data relating to 
other magnetic effects* 

20 Preferably the trend of the gradient of the ship bias is 

determined by applying the raw magnetic gradient data to a 
B-Spline filter with smoothness 0.8. 

Preferably the processing comprises the application of the 
25 corrected gradient data over sample intervals, integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low pass filter to 
the integrated total magnetic intensity data to provide 
the data output. 

30 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
the following equation: 

35 TMI= 2 [GcAxiiM+Mw 

i=i 

wherein G c is the corrected gradient data 
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obtained after subtraction of the trend of the gradient of 
the ship bias from the raw gradient data, and 

M st at is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other . 



Preferably the gradient of the ship bias is calculated 
according to the following equation: 

10 

Gt>ias = M f^~ M ^ = IDj^-DWHMM-MM] 
A/ Al 

where : 

15 M/x) = M e (x) + £>(/,) + M b (ti) 

MApc) = M e (x) + D(t 2 ) + M b (t 2 ) 



20 



where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
MoOc), diurnal variation D(t x ), ship bias tMti) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time t L and along track distance x, 
and at some later time, t 2 , D(t 2 ) is the diurnal variations 
25 sensed by the rear sensor, and 

Mb(t 2 ) is the ship bias at time t 2 detected by the 
rear sensor. 



30 



The invention may also be said to reside in a method of 
obtaining gradient data for an exploration program, the 
method including: 

towing first and second sensors behind a ship 
along predetermined survey lines; 

obtaining raw magnetic gradient data from the 

35 sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
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subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a data output. 

Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

Preferably the trend of the gradient of the ship bias is 
determined by applying the raw magnetic gradient data to a 
B-Spline filter with smoothness 0.8. 

Preferably the further processing comprises the 
application of the corrected gradient data over sample 
intervals, integration of the corrected gradient data into 
total magnetic intensity data, and application of a low 
passed filter to the integrated total magnetic intensity 
data to provide the data output. 



Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
25 the following equation: 

TMI=£ [G c A^,(0]+M ftof 
i=i 

wherein G c is the corrected gradient data 
obtained after subtraction of the ship bias trend from the 
raw gradient data, and 

Motat is the total magnetic intensity value at a 
point .representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other. 

Preferably the gradient of the ship bias is calculated 

H.NLuisaVKeepXSpeciXBHP Prov - Method ol Proceoain, Marine Miotic Orient Decdoc 3/12/01 
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according to the following equation: 

= M f {x)-M r {x) _ [P(r, ) - Dih )] + [M, (f, ) = (r 2 )] 
A/ A/ 

5 where : 

M/*) = M e (x) + D{t x ) + M*(f 0 
M^x) = M*(*) + Z>(* 2 ) + Aft(r 2 ) 

10 

where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(ti), ship bias M b (ti) caused by 
ship-induced field, sensor drifting and sensor 
15 perturbation, all at time ti and along track distance x, 

and at some later time, t 2 , D(t 2 ) is the diurnal variations 
sensed by the rear sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
rear sensor. 

20 

The invention may also be said to reside in a method of 
subsurface exploration to determine viability of drilling 
in a marine, environment by considering magnetic data 
relating to the environment, and wherein the magnetic data 
25 has been obtained by: 

obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
30 subtracting the trend from the raw magnetic 

gradient data to obtain corrected gradient data;. and 

processing the corrected gradient data to provide 
a data output. 

35 Preferably obtaining the magnetic data includes the step 
of determining the gradient of the ship bias from data 
obtained by the sensors, and determining the trend of the 
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gradient of the ship bias from the gradient of the ship 
bias . 

Preferably the trend of the gradient of the ship bias is 
determined by applying the raw magnetic gradient data to a 
B-Spline filter with smoothness 0*8. 

Preferably the processing comprises the application of the 
corrected gradient data over sample intervals, integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low passed filter to 
the integrated total magnetic intensity data to provide 
the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
the following equation: 



20 wherein G c is the corrected gradient data 

obtained after subtraction of the ship bias trend from the 

raw gradient data, and 

Mfitat is the total magnetic intensity value at a 

point representing the start of the survey or at a 
25 location where a survey line and a tie line cross each 

other • 

Preferably the gradient of the ship bias is calculated 
according to the following equation: 



Gbias = M /(*>~^(*) ^ [D(t { )-D(t 2 )} + [M b (t { )-M b (t 2 )} 
&l Al 



where : 

3 5 Mj(x) = M e (x) + D(t t ) + M b {t x ) 
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M&c) = M e (x) + D(t 2 ) + M b (t 2 ) 

where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
5 Me(x), diurnal variation D(ti), ship bias M b (ti) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time t x and along track distance x, 
and at some later time, t 2 , D(t 2 ) is the diurnal variations 
sensed by the rear sensor, and 
10 M b (t 2 ) is the ship bias at time t 2 detected by the 

rear sensor. 



The invention may still further be said to reside in a 
method of drilling for a deposit in a marine environment, 
1 5 inc luding ; 

determining the location of drilling from data 
which has been obtained and which indicates the possible 
existence of the deposit; and 

which location is also determined by magnetic 
20 data which has been obtained by towing magnetic sensors 
behind a ship, the magnetic data being processed by: 

obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
25 bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; 
and 

processing the corrected gradient data to provide 
30 a data output* 



Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
3 5 determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 
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Preferably the trend of tlie gradient of the ship bias is 
determined by applying the raw magnetic gradient data to a 
B-Spline filter with smoothness 0.8. 

5 Preferably the processing of the corrected gradient data 
comprises the application of the corrected gradient data 
over sample intervals, integration of the corrected 
gradient data into total magnetic intensity data, and 
application of a low passed filter to the integrated total 
10 magnetic intensity data to provide the data output. 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
the following equation: 



wherein G c is the corrected gradient data 
obtained after subtraction of the ship bias trend from the 
raw gradient data, and 



point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other . 

25 Preferably the gradient of the ship bias is calculated 
according to the following equation: . 



15 



stat 



20 



M s tat is the total magnetic intensity value at a 



Gbias = 




30 



where : 



M/x) = M e {x) + D{t x ) + M h (t0 



MM = M e (x) + D(t 2 ) + M b (t 2 ) 



35 



where M f is the magnetic field as measured by the 
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front sensor and which consists of the environmental field 
Me(x), diurnal variation D(ti), ship bias M b (t x ) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time t x and along track distance x, 
5 and at some later time, t 2 , D(t 2 ) is the diurnal variations 
sensed by the rear sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
rear sensor. 

10 Brief Description of the Drawings 

A preferred embodiment of the invention will be described, 
by way of example, with reference to the accompanying 
drawings in which: 

Figure 1 is a view showing a marine magnetic data 

15 survey; 

Figure 2 is a flow chart showing the preferred 
method of processing the marine magnetic gradient data 
according to the preferred embodiment of the invention; 

Figure 3A and Figure 3B are graphs showing data 

2 0 obtained according to the survey and as processed 

according to the preferred embodiment of the invention; 

Figure 4 is a representation of a survey using 
. data created according to a prior art technique; and 

Figure 5 is a diagram similar to Figure 4, but 
25 using data obtained according to the preferred embodiment 
of the present invention. 

Description of the Preferred Embodiment 

With reference to Figure 1, a marine survey to acquire 

3 0 magnetic data is illustrated. A ship 10 tows first 

magnetometer Ml and second magnetometer M2 by a tow line 
12. The ship travels along survey line 14 and the line 
16, which are dotted, represent the ship-induced magnetic 
field. 

35 

In Figure 1, it is assumed the ship is in the southern 
hemisphere and is heading north-northwest, and the sensors 
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are drifting towards tlie northeast. In this case, the 
recorded gradient data between sensors Ml and M2 will have 
a larger distortion than if the sensors were drifting 
towards the southwest . 

5 

Ideal conditions would assume that the sensors Ml and M2 
are towed directly behind the ship on the survey line and 
that therefore, after a given time period, the sensor M2 
would occupy the same position as was previously occupied 

10 by the sensor Ml. However, due to drifting of the sensors 
because of ocean currents and the like, this ideal 
situation in general cannot be achieved. Nevertheless, 
the processing methods used to process magnetic data do 
assume that the sensor M2 will occupy the same position as 

15 the sensor Ml, previously occupied, but at a later time. 
The processing method according to the preferred 
embodiment enables ship bias to be accounted for much more 
accurately than any prior art techniques and therefore 
enables the sensors to be towed closer to the ship and to 

20 be separated by a smaller distance. This therefore 

provides much greater accuracy in the measurements because 
the amount of drift is not as great due to the shorter tow 
line and also because the separation between the sensors 
is smaller, the likelihood of the rear sensor occupying 

25 the same position as the front sensor at a later period in 
time is greatly increased. 

At any time t x and along-track distance x, the front sensor 
measures magnetic field - Ajr consisting of the environment 
30 field Me(x), the diurnal variation D(t x ) , ship bias Jfe(fci) 
caused by ship induced field, sensor drifting and sensor 
perturbation. 

Mj(x) = M e {x) + D(t ,) + M b (t { ) 

35 

At some later time t 2 , the rear sensor makes a measurement 
at the same spatial point: 
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MAx) = M e (x) + D(t 2 ) + M b (t 2 ) 

The difference between the two measurements at the same 
5 location is: 

Mj{x) - M^x) = [D(t x ) - D(t 2 )] + [M b {h) - M b (t 2 )] 
Ship bias gradient (Step 2, Figure 2) can be expressed as: 

10 

_ M f (x)-M r (x) _ [D(f, ) - £>(/ 2 )] + [M h (t l )-M h (f , )] 

^ bias — — — — — 

A/ Al 

It should be mentioned that the calculated ship bias 
gradient is still affected by diurnal variation. In 

15 practice we notice that the sea current drifting is a kind 
of long wave-length variation so ship bias should be 
coincident with such current drifting. Step 2 in Figure 
2, which is the calculation of the gradient of the ship 
bias # is actually an estimate of the gradient of the ship 

20 bias which is obtained from the raw data collected by the 
sensors Ml and M2 . The data collected by the sensors will 
include raw magnetic gradient data which contains many 
signals including environmental magnetic signals, diurnal 
signals and the ship bias, as well as instrument bias and 

25 drift. The trend of the gradient. of the ship bias is 

obtained from this estimate by applying the estimate of 
the gradient of the ship bias to a B-Spline filter, as 
will be described in more detail hereinafter. Therefore 
the trend of the gradient of the ship bias (Gtrend) (Step 

3 0 3, Figure 2) is used for the correction along a track 
line. 



The gradient correction of the ship bias (GO can be 
expressed as follows: 

35 

G C = G- G tren d (Step 4, Figure 2) 
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As said before such drifting can cause large bias effects. 
Consequently, when the sensor drifting and perturbation 
take place simultaneously the bias zig-zag will be around 
5 its mean values. The variable mean value of the bias 
along a survey line is regarded as bias trend since the 
high frequencies of sensor perturbation only create random 
noise around the bias trend and after integration its 
effect can be eliminated from raw gradient data. 

10 

After removal of ship bias effects TMI (Total Magnetic 
Intensity) can be calculated (Step 6) through integration 
of magnetic gradient data: 

I 

15 TMI= « [GcAxi(t)] + Mstat 

Here, AsTi(t) is the sampling distance along the survey 
line. Matat is TMI value at a point of survey start or at 
the location where the survey line and a tie line cross 
20 each other. 

. Figures 3A and 3B represent an actual example of the 
preferred embodiment of the invention which relates to raw 
data obtained from a known region. 

25 

In Figure 3A, the trace 20 represents the total ship bias. 
The line 21 represents the bias trend and it can be seen 
from the left hand side of the trace 20 to the right hand 
side of the trace 20, the bias trend 21 changes non- 
30 uniformly about the value of -0.08 on the graph of Figure 
3A. A particular value for the bias trend can therefore 
be obtained for various intervals and this value can be 
subtracted from the raw data in order to provide the 
corrected data. The trend of the gradient of the ship 
3 5 bias is preferably determined by calculating the gradient 
. of the bias data 20 and then applying that gradient data 
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to a B-Spline filter with smoothness 0.8 in order to 
provide the representation of the ship bias 21. 

In the preferred embodiment of the invention, the gradient 
5 of the bias data 20 is calculated simply by using the raw 
magnetic gradient data obtained from the sensors because 
when that data is applied to the B-Spline f ilter, all that 
remains is the gradient of the ship bias component. The 
ship bias 21 can then be subtracted from the calculated 
10 raw gradient. 

The trend of the gradient of the - ship bias is a non-linear 
function which is represented by the line 21 in Figure 3A. 
The line is a measure of how the bias changes over time as 

15 the sensors are towed behind the ship. As is clearly seen 
from line 21 , the trend is not constant or merely an 
average, but rather represents the fluctuation of the bias 
gradient and at some times, is above the value of 0.08 and 
sometimes below that value. As noted above, the trend is 

20 determined by applying the gradient of the ship bias to a 
B-Spline filter with smoothness 0.8. However, in other 
embodiments, the filter could have different smoothness 
depending on the region from which data is collected and 
the nature of the data which is collected. In general, 

25 the function of the filter is to smooth the trace 20 so 

that some meaningful value of the bias trend at particular 
time periods can be obtained. Thus, the filter is 
effectively determining the peaks and troughs of the trace 
20 and fitting a curved line between those peaks and 

3 0 troughs, which gives a measure of how the gradient of the 
ship bias is changing with time as the sensors are towed 
behind the ship. 

In the example of data obtained from the known region, 
3 5 dashed line 30 represents total magnetic intensity data 
which is provided according to the prior art processing 
technique. The line 32 represents the integrated total 
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magnetic intensity data from the bias corrected gradient 
data according to the preferred embodiment of the 
invention. The line 34 represents the diurnal variation 
at a station some 500 km away from the survey area and 
5 line 36 is the observed field data including the diurnal 
effect . 

It can be seen from the conventional processing technicjue, 
which gives line 30, false anomalies are provided which do 
10 not show up on the line 32 , which is produced according to 
the present invention. 

In general, the results of the magnetic survey are 
produced in a colour chart. Black and white 
15 representations of a colour chart are shown as Figure 4 in 
relation to conventional processing of data obtained from 
a known region and Figure 5 by a processing method 
according to the preferred embodiment of this invention. 

20 The characteristics of the Igneous basement of the known 
region are well known and it can be seen that data 
processed according to the present invention produces a 
much clearer indication of actual magnetic structures than 
Prior art technique which include considerable 

25 interference and data which may mislead analysers into 

considering that magnetic structures exist which are not 
actually there, or that magnetic structures do not exist 
when a magnetic structure actually is present. 

3 0 Thus, the preferred embodiment of the present invention 
provides data which more accurately reflects the 
likelihood of magnetic structures which can then be used 
as a basis for determining the viability of a drilling 
operation in an exploration or recovery program. 

35 

Since modifications within the spirit and scope of the 
invention may readily be effected by persons skilled 
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within the art, it is to be understood that this invention 
is not limited to the particular embodiment described by- 
way of example hereinabove. 
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